Structure-dependent exchange in the organic magnets Cu(II)Pc and Mn(II)Pc 
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We study exchange couplings in the organic magnets copper(II) phthalocyanine (Cu(II)Pc) and 
manganese(II) phthalocyanine (Mn(II)Pc) by a combination of Green's function perturbation the- 
ory and ab initio density-functional theory (DFT). Based on the indirect exchange model our 
perturbation-theory calculation of Cu(II)Pc qualitatively agrees with the experimental observa- 
tions. DFT calculations performed on Cu(II)Pc dimer show a very good quantitative agreement 
with exchange couplings that our theoretical group extracts by using a global fitting for the mag- 
netization measurements to a spin-| Bonner-Fisher model. These two methods give us remarkably 
consistent trends for the exchange couplings in Cu(If)Pc when changing the stacking angles. The sit- 
uation is more complex for Mn(II)Pc owing to the competition between super-exchange and indirect 
exchange. 

PACS numbers: 71.10.-w,71.15.Mb,71.35.Gg,71.70.Gm,75.10.Pq,75.50.Xx 



I. INTRODUCTION 

Recently molecular spintronics [l|, 0, d, 0, [f| has be- 
come a very active inter-disciplinary topic. This is be- 
cause localized spins in molecular complexes can have 
very long spin relaxation times (up to of order one sec- 
ond) @, while the chemical engineering of such com- 
plexes is much more flexible than is the case in con- 
ventional inorganic-semiconductor electronics. The main 
building blocks of molecular spintronics, namely radicals 
containing localized electrons, are promising candidates 
both for spintronics and for quantum information pro- 
cessing. Against this background, experimental and the- 
oretical studies of the magnetism in spintronics-related 
organic materials are crucial for the development of de- 
vices such as molecular magnetic random-access memory 
(MRAM) [1,0,1. 

There is a long history of research on metal phthalo- 
cyanines (MPc) , because of their commercial applications 
and excellent electro-optical properties @. In particular, 
their magnetic properties have been extensively studied 
[IS M, QL Mn(II)Pc was one of the first molecular 
magnets its properties were shown to depend criti- 
cally on the stacking of the planar molecular 7r-systems. 
We define a stacking angle as shown in Figure [TJ the (3- 
Mn(II)Pc crystal (stacking angle 45°) was found to be 
ferromagnetic, while the a-Mn(II)Pc thin film (stacking 
angle 65°) was shown to be antiferromagnetic. However, 
there is still debate about the details of the molecular 
stacking in the a-phase Cu(II)Pc material: transmission- 
electron diffraction (TED) observations for Cu(II)Pc on 
a KC1 (001) surface suggested that the stacking ori- 
entation in the a-phase was the so-called " x " model as 




FIG. 1: Schematic arrangements of (a) the a- and (b) /3-phase 
MPc structures. The stacking angle is (j> in each case. 



shown in Figure [2ja). However, the most recent TED 
experiments |l5j| by contrast indicated the orientation is 
the " +" model shown in Figure [D(b) . The arrows in 
Figure [5] show the direction of the displacement of the 
Pc molecules between successive layers; these directions 
differ by a rotation of 45°. Since this issue is not yet 
resolved, and given that the substrate used in [TJ, [l5| 
differs from that used in [l6l |. in the following discussion 
we adopt the "+" model. 

Recently Heutz [l6j], et ai., have performed further 
magnetic measurements on different phases of Cu(II)Pc 
and Mn(II)Pc by using SQUID magnetometry. In the 
remainder of this paper, we describe these spin systems 
by using a Heisenberg spin-chain model [13] , which is be- 
lieved to be a good description for organic systems con- 
taining localized spin centers: 
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H e ff = — 2J Si ■ Sj. 



(1) 



Note that with the sign convention we adopt, a positive 
exchange constant J corresponds to ferromagnetic cou- 
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(a) 




(b) 



stants. 

These experiments clearly show the switching of mag- 
nitudes and signs of the exchange couplings as the molec- 
ular packing varies from phase to phase, and also that the 
magnetic properties are determined by the structure (a 
versus /?), not by the sample morphology (powder ver- 
sus thin film). The results confirm the previously mea- 
sured [l(| E3 difference between the a and (3 phases of 
Mn(II)Pc, and also show that a corresponding difference 
exists for Cu(II)Pc, although in this case the (3 phase is 
paramagnetic rather than ferromagnetic. 

Despite the long history of experimental work on MPc, 
there have been very few systematic theoretical studies of 
the mechanisms underlying the variation in the exchange 
interactions; the problem is complicated by the molecu- 
lar structure and rather weak spin-spin interactions com- 
pared to conventional inorganic semiconductors. In this 
paper we aim to gain both a picture of the physics driving 
the structure-dependent exchange, and a quantitative un- 
derstanding of its magnitude, in Cu(II)Pc and Mn(II)Pc. 

Our remaining discussion falls into five sections. In 
jfiTl we introduce the different mechanisms for exchange 
and describe state-of-art quantitative methods to eval- 
uate exchange interactions by using DFT and the bro- 
ken symmetry concept. We also describe the atomic and 
electronic structure of the systems we consider. In 
we perform Green's-function perturbation theory calcu- 
lations for exchange interactions to get a rough picture 
of essential physics. In EjlVI we use DFT with the hy- 
brid exchange-correlation functional B3LYP to evaluate 
the exchange interactions quantitatively. At the end we 
draw our conclusions in SjVl 



FIG. 2: Two models for stacking orientations for a-Cu(II)Pc: 
(a) "x"; (b)"+". 



II. THEORETICAL OVERVIEW 



pling, while a negative J describes anti-ferromagnetic in- 
teractions. In these experiments Mn(II)Pc powder sam- 
ples (/3-phase) show strongly ferromagnetic (FM) cou- 
pling with J 11.45K, while Mn(II)Pc films (a-phase) 
grown on an inert Kapton substrate shows a relatively 
much weaker anti-ferromagnetic (AFM) coupling with 
J w — 1.61K. Cu(II)Pc powder (/3-phase) is found 
to be very weakly ferromagnetic (indeed nearly para- 
magnetic) with J OK, but Cu(II)Pc films (a-phase) 
and Cu(II)Pc whose growth is templated by a layer of 
3,4,9,10-perylenetetracarboxylic dianhydride (PTCDA) 
pre-deposited on the Kapton substrate are found to be 
more strongly anti-ferromagnetic (J sa — 1.50K). The 
exchange constants for Mn(II)Pc are extracted from the 
intercept of the inverse susceptibility versus temperature, 
while those for Cu(II)Pc are found by a global fit of the 
experimental data to a finite S = 1/2 Heisenberg spin 
chain model (the so-called Bonner-Fisher model [la]). 
This model is expected to be sufficiently accurate, de- 
spite its neglect of inter-chain couplings, provided the 
temperature is not too low relative to the exchange con- 



A. Qualitative description of exchange interactions 
between localized electrons 

The various processes contributing to the exchange in- 
teractions between localized spins were extensively con- 
sidered by Anderson (l9| . His original paper consid- 
ered exchange interactions between magnetic ions in ionic 
crystals, but the same concepts and arguments apply in 
the case of the bi-radical studied here. The direct ex- 
change interaction originates from a quantum exchange 
term of the Coulomb interaction between localized elec- 
trons, e.g., d-electrons in a magnetic ion; this always gives 
rise to ferromagnetic exchange interactions. However, for 
MPc the direct exchange interaction can generally be ne- 
glected owing to the large molecular size and the local- 
ization of the metal electrons. 

The swper-exchange arises from the terms in the 
Hamiltonian that tend to delocalize electrons. It is then 
necessary to take the hopping perturbation to higher or- 
der in order to reach an excited state in which an elec- 
tron is transferred onto a neighboring magnetic site. For 
example, in a simple Hubbard model the second order 
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in perturbation theory leads to an exchange interaction 
—2t 2 /U, where t is the transfer integral between sites and 
U is the on-site Coulomb repulsion. In a more realis- 
tic model there are many more super-exchange pathways 
but the same principles remain: according to the defi- 
nition we adopt, for a process to be classified as super- 
exchange it should proceed through virtual states involv- 
ing the migration of charge between the magnetic centers. 
The super-exchange mechanism is especially important 
when the magnetic atoms are separated by non-magnetic 
species, for example inorganic anions or organic ligands, 
through which charge can pass from one magnetic atom 
to another. We should note that the super-exchange van- 
ishes not only when the distances are large but also when 
the the transfer of electrons between the magnetic centers 
and these covalent "bridges" is symmetry- forbidden. 

The indirect exchange interaction between electronic 
spins is similar to the indirect exchange interaction be- 
tween nuclear moments, which is mediated by the con- 
duction (or valence) electrons: a polarization of the con- 
duction electrons around one local moment is propagated 
to another, giving rise to an effective interaction. This 
nuclear interaction was first discovered by Ruderman and 
Kittel [2lj . and independently in molecular physics by 
Ramsey and Purcell [13, and by Bloembergen and 
Rowland [24J; the generalization to localized electronic 
moments is due to Kasuya [25[ and Yosida [26| . In metals 
it leads to a long-range exchange coupling that is oscilla- 
tory in sign. 

We can sharpen the distinction between the different 
types of exchange by writing the full electronic Hamilto- 
nian as 

H H^pin ~l~ -^ligand ~l~ ^direct ~l~ ^hop ~t~ ^polarize- (2) 

The terms are defined as follows. H spin corresponds to 
the isolated spins, and H\i gan d to the rest of the material 
(decoupled from the spins), both being taken to include 
an appropriate (spin-independent) mean field. We take 
Hq = -ff S pm + -ffiigand to be our unperturbed Hamiltonian. 
Vdi rec t then involves the direct Coulomb interaction, cou- 
pling the spins to one another without involving the lig- 
ands; Vh p involves all processes that transfer an elec- 
tron between the magnetic species and the ligands; while 
polarize includes all other (non-charge-transferring) in- 
teractions between the magnetic species and the ligands. 

We can then develop a perturbation expansion for the 
full Green's function G (as, for example, in 20]) as 

G = G + GqVG + GqVGqVGo + ... (3) 

where V is the perturbation Vdircct + Vh op + ^polarize and 
Go is the Green's function corresponding to unperturbed 
Hamiltonian. The effective exchange is then recovered 
by computing an effective Hamiltonian 

E-H cS = G- 1 (4) 

within a ground-state manifold of configurations differ- 
ing only in the spin orientations. Within this picture, 



the direct exchange corresponds to the first-order term 
involving Vdircct, while super-exchange and indirect ex- 
change correspond to higher-order terms involving Vhop 
and V^oiarizo respectively. It is clear from the definitions 
of the various terms in V that the virtual states that cou- 
ple to the ground-state manifold are orthogonal; there- 
fore (at least to low orders in V) we do not need to con- 
sider cross-terms between the different operators. 

In this paper we will neglect Vdircct, for the rea- 
sons given above. Appropriate expressions for Vhop and 
Vp i ar i Z o are given in EjTTT] below. 

B. Quantitative calculation: DFT calculations of 
exchange couplings in bi-radicals 

Density Functional Theory (DFT) [H [H, [H M, 
l32l HH is a powerful tool for accurate prediction of the ex- 
change interactions in chemically complex wide-gap ma- 
terials. However current density functionals, based on 
Kohn-Sham theory, give a poor representation of sin- 
glet states containing a pair of localized electron spins 
since the Kohn-Sham orbitals are constrained to re- 
spect the symmetry of the system and are therefore 
generally formed from linear combinations of the single- 
center wavefunctions. One has to use instead a so-called 
'broken-symmetry' method 29], in which the magnetic 
orbitals are localized in different radical centers, with 
their spins oppositely-aligned. Recently, Martin and Illas 
[32l checked the performance of different exchange- 
correlation functionals in the calculation of magnetic cou- 
plings, and found the choice of exchange functionals is 
extremely important, while the role of the correlation 
functional is minor. Although the precise reasons are un- 
clear, it is empirically found that it is necessary to mix 
some proportion of exact exchange into the functional 
in order to obtain results that agree with experiment or 
with higher-quality quantum chemistry results for small 
molecules — crudely, this may be understood as requiring 
some balance between the over-localization of electrons 
in a Hartree-Fock calculation, and the excessive dereal- 
ization in standard density functionals. In particular the 
B3LYP functional [H, Ha, l37| , which mixes about one 
quarter Hartree-Fock exchange, has been found to give 
good results for di-nuclear molecules, organic bi-radicals, 
and spins localized at defects in carbon-containing ma- 
terials HI, HI Hi- In SB we perform DFT with 
B3LYP exchange-correlation functional to calculate ex- 
change interactions in Cu(II)Pc dimers based on this 
broken-symmetry concept. 

Although the DFT and perturbative approaches to the 
problem appear at first sight to be quite different, one can 
think of them as representing in different ways the same 
response of the electronic system to its spin-dependent 
interactions. In the case of the DFT, this response is 
represented as a change in the Kohn-Sham states, while 
in the model approaches the many-electron wave function 
responds by including small admixtures of excited-state 
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configurations. 



C. The electronic structure of isolated Cu(II)Pc 
and Mn(II)Pc 

Before we perform the perturbation theory calcula- 
tion for the exchange interaction, we need to under- 
stand the nature of the one-electron states in the isolated 
molecules. We used the Gaussian 98 code (38|, perform- 
ing a DFT calculation with the B3LYP [37] exchange- 
correlation functional and a 6-31G [3^| basis set, to op- 
timize the molecular geometry of isolated Cu(II)Pc and 
Mn(II)Pc molecules. We then use the key Kohn-Sham 
states emerging from DFT calculation which are near- 
est to the HOMO-LUMO gap as a basis to perform a 
separate perturbation-theory calculation. 

Our level scheme for Cu(II)Pc is shown in Figure (3^a). 
The states are labeled by the irreducible representations 
of the D^h point group. From the Mulliken population 
analysis we can identify b\ g as a metal d-orbital which is 
hybridized with the Pc ring (Mulliken charge 0.30). The 
total Mulliken charge on copper is +0.97. The total Mul- 
liken spin density on the copper atom is 0.68; this is con- 
sistent with the existence of one singly-occupied orbital, 
with a spin mainly but not entirely localized on the cop- 
per atom. The overall symmetry of the electronic state 
is 2 Big. We found that the occupied molecular orbital 
with the largest Kohn-Sham eigenvalue is not the singly- 
occupied big state, but the a\ u state; this is slightly dif- 
ferent from the early extended Huckel calculations [io| . 
This highlights the importance of two-electron Coulomb 
terms in determining the configuration: doubly occupy- 
ing the big state would incur a large Coulomb penalty 
because the charges would spend much of their time local- 
ized in the Cu 3d states, whereas the double-occupancy 
penalty for the more diffuse ai u state is much smaller. 

For Mn(II)Pc, the expected total spin is S = 3/2 
[ToL Hf3 |. Our Gaussian calculation gave the overall 
electronic configuration 4 Ai g , with three singly-occupied 
one-electron levels having ai g and e g (twice) symmetry. 
The total Mulliken charge on Mn is +1.14; again this is of 
the same order as, but somewhat less than, the nominal 
+2 valence. The Mulliken spin density on the Mn atom 
is +3.1. Our calculation results are similar to those in 
which give ^Ai g total symmetry and the same single- 
occupied states. However, this picture of the electronic 
structure is not the only one. Liao et. al [4l[ also used 
DFT methods and found an electronic configuration in 
which the three unpaired electrons occupy the ai g , &2g, 
and one e g state while the other e g state is doubly oc- 
cupied, to give an electronic symmetry 4 E g . This cal- 
culation is in agreement with the more recent magnetic 
circular dichroism (MCD) and UV-vis measurements [42j 
of the molecule in an argon matrix, but differs from the 
early magnetic measurements of solid Mn(II)Pc. It is 
possible that the 4 E g configuration (which would lead to 
a Jahn- Teller distortion, because of its orbital degener- 
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FIG. 3: The schematic of the key states from our Gaus- 
sian DFT calculations of the electronic structures of isolated 
(a) Cu(II)Pc and (b) Mn(II)Pc. The majority-spin contri- 
butions to the Mulliken charge on the transition-metal atoms 
are shown in brackets (the minority-spin contributions for the 
doubly-occupied states are similar). 

acy) may be favored in the isolated molecule or in the 
argon matrix, with the 4 A2 g state favored in the bulk 
material (where no Jahn- Teller distortion has been ob- 
served). 

III. GREEN'S FUNCTION PERTURBATION 
CALCULATIONS 

A. Super-exchange calculation 

We aim to understand the mechanism of exchange 
couplings between neighboring Mn(II)Pc and Cu(II)Pc 



molecules observed in experiments [HI Ell • We first 
consider the super-exchange contributions. 



1. Cu(II)Pc 

As explained in CTfl the super-exchange contribution 
is generally dominant when considering the exchange in- 
teraction between localized spins in insulating materials. 
However, Cu(II)Pc is an example of a situation where 
this interaction is expected to be negligible. This is be- 
cause the unpaired spin is located in a bi g orbital [io| . 
but there is no low-energy state of b\ g symmetry in the 
ligand available to hybridize with it. Therefore, as long 
as the symmetry of the molecule remains D^h (believed 
to be an excellent approximation even in the crystal) the 
spin-carrying electron is "tied" to the Cu site and no 
super-exchange can take place except by direct hopping 
from the Cu orbitals onto the neighboring molecule. The 
amplitude for this process is expected to be very small. 



2. Mn(II)Pc 

In the case of Mn(II)Pc, there are three unpaired 
electrons per molecule occupying the a\ g , e gx , and e gy 
molecular orbitals (see £111 C|) . We can call these metal 
states because these orbitals originate from the splitting 
of atomic 3d-states in the molecular environment which 
has Dih symmetry. The two S — 3/2 spins of the indi- 
vidual molecules can be combined to form a total spin of 
0, 1, 2, 3. We therefore need, in principle, three indepen- 
dent parameters in the spin Hamiltonian to characterize 
fully the relative energies of these states. If we neglect 
spin-orbit coupling, the Hamiltonian must be invariant 
under simultaneous rotations of both spins and therefore 
must have the form 

H = -Ji(S A ■ S fl ) - J 2 (S A ■ S B ) 2 - J 3 (S A ■ S B ) 3 , (5) 

where J\ , J2, and J3 are exchange couplings and A, B 
label the molecules. We can find all three parameters 
from the 4x4 Hamiltonian matrix spanned by the 
states with total z spin angular momentum Ms = 0: 
|(3/2U, (-3/2) B ), |(1/2U, (-l/2) fl ), K-1/2U, (l/2) s ), 
and |(— 3/2)a, (3/2)b). We use a similar method 
to that described in [2Cj : we construct the effective 
Hamiltonian based on an extended Hubbard model by 
Green's-function perturbation theory, compare this with 
equation ([5]), and extract the exchange constants. For 
simplicity we include only intermediate states where a 
single electron is transferred between adjacent Mn ions 
via the ligand e g states (i.e. we neglect the possibility 
that two or more electrons transfer together). We also 
neglect direct electron transfer between the d states 
of the Mn ions, because these states are quite well 



localized. Our extended Hubbard model reads: 
H = H + H t + H p : 

h q = EjUj 

i 

+V pmn (nE gx A + n Egy A)(n egac A + n CgyA + n aigA ) 
+u gx h aig A{h egx A + n egyA ) 
+u xx (n egx Aine gac Ai + n egvA ^n egy Ai) 
+u gg n aig Afn aig Ai} 
+{A^B}; 

Ht = {^2[ t prn a cl aigA)a (c {Eg:c B)a+C { E gv B)a) 
a 

+tp m xl{c\ egiA)a C { E^B)a + c\ egyA)a C (EgyB )a) 
+ t pmx2(c\ egxA)lT C { E gv B) < j + c[e g „ A)a d (E g , B)a) + h.C 

+{A^B}; 
H p = -Ks ■ S. 

Here ai g ,e gx , and e gy label the metal states; E gx and 
E gy label LUMO ligand states for distinguishing metal 
and ligand e g states [9]. Hq includes the single-particle 
energies Ei where i runs through a% g , e gx , e gy , E gx , E gy , 
the Coulomb interaction between metal and ligand states, 
and the on-site Coulomb interactions. There are two 
parts in the perturbation: one is H t which transfers elec- 
trons between molecules and the other is H p represent- 
ing the interaction between i-spin (s) in the ligand and 
|-spin (S) on the metal within the molecule. We sup- 
pose that Hp is itself ultimately a representation of a 
further super-exchange processes, and therefore like H t 
originates in Vhop as defined in §11 Al t pmg , t pmx i and 
t P mx2 are the inter-molecular hopping integrals shown in 
Fig. [4] Eg and e g are the energies of ligand and metal 
e g states relative to the energy level of a\ g state, u gx is 
Coulomb interaction between the a\ g and e g levels, u xx 
is the Coulomb interaction between two degenerate Mn 
e g states, and v pmn is the Coulomb interaction between 
the Mn and Pc states within a molecule. 

From this Hamiltonian we can see when one elec- 
tron is transferred from the metal state of molecule 
A to the ring state of molecule B where it can in- 
teract with the Mn spin; it is through the interaction 
Hp that the spin projections ttla and raj associated 
with the two molecules can change, thereby coupling the 
four spin states: |(3/2) A , (-3/2) B ), |(1/2) A , (-1/2) B ), 
|(-1/2) A ,(1/2) B ), and |(-3/2) A , (3/2) B ). Note that in 
D^h symmetry, the e g states of Mn can hybridize ef- 
fectively with the e g states of the ring; the existence 
of unpaired spins in the e g states is what makes super- 
exchange processes much more important in the case of 
Mn(II)Pc. 

There are 25 spatial configurations and each has four 
possible spin states, giving a total of 100 states. We 
construct the 100 x 100 Hamiltonian matrix for Hq and 
V, and then extract the effective Hamiltonian within the 
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4x4 low-energy subspace by using Green's function per- 
turbation theory [2(| to calculate the energy shifts. By 
comparing this low-energy subspace with the |-spin cou- 
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FIG. 4: The possible inter-molecular transitions. 

We note several features of this result. First, the dom- 
inant terms are those proportional to i 2 , i.e. where elec- 
trons are exchanged once between the molecules, as ex- 
pected in a super-exchange process. Second, the leading 
term in Ji is proportional to Kt 2 , in J2 to K 2 t 2 , and 
in J3 to K 3 t 2 ; this is because H p only couples states 
in which rriA and ms alter by one unit of angular mo- 
mentum. Finally, assuming the Coulomb energies are 
all large and positive, J\ is always the same sign as K, 
irrespective of the values of the various hopping terms. 
In general we expect that K, since it is dominated by 
super-exchange, will be negative (corresponding to anti- 
ferromagnetic coupling in our sign convention) and there- 



pling matrix, we find that it can be written in the form 
(p)]) . with parameters 



(10) 



(11) 



(12) 

I 

fore J2 will also lead to anti-ferromagnetic coupling in- 
dependent of the orientation of the molecules. 

Our conclusion about the failure of the super-exchange 
interaction to change sign contrasts sharply with the ex- 
planation given by Barraclough et al. [10| and by Ya- 
mada et al. (llj for their experimental results, which 
they ascribe to the competition between different super- 
exchange pathways operating via nitrogen atoms. How- 
ever, this argument fails to take into account correctly 
the spin algebra — in particular, it ignores the fact that 
the three electron spins on each Mn atom are in fact tied 
together via strong intra-atomic Coulomb interactions, 
and so cannot be flipped independently. 



B. Indirect exchange calculation 

1. Gu(II)Pc 

For the indirect exchange scheme in the Cu(II)Pc 
dimer (Figure O , the unpaired electron spin on the metal 
polarizes the ligand by the two-body Coulomb interac- 
tion; this spin polarization can transfer to the neighbor- 
ing molecule by orbital hybridization, and there inter- 
acts with the unpaired spin of the neighboring molecule's 
metal ion. 

Because the LUMOs are e g ligand states, we should 
consider the filled ligand states with the same symmetry. 
The two-body Coulomb interaction can be represented in 
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metal spin operators and ligand spin-polarization opera- 
tors owing to the preservation of total S z in the isolated 
molecule. We label the spatial LUMO state of the ligand 
by using "X", the filled states "G", and metal b\ g state 
" b" . We use the following notation for the two-electron 
integrals: 



{a,b\c,d}= [ drdr'a(r)*b(r')* , 1 . cfrW). (14) 
J |r-r'| 

We can now apply Green's-function perturbation the- 
ory [2(3] to this problem; the perturbation includes the 
Coulomb interaction v which can polarize the spin in 
a ligand, and hopping t that transfers this polarization 
from one molecule to another. We find that the leading 
term J\ in the spin-i couplings is given by: 



FIG. 5: Cu(II)Pc electron configuration and indirect exchange 
scheme diagram, 
two empty e g states(LUMO). 



second-quantized form as 

nn' A R n K ' I 



aa' ABDE 



(13) 



where c and c< are the electron annihilation and creation 
operators, and A, B, D, E may each represent a metal or- 
bital or a ligand orbital. Since we wish to consider pro- 
cesses in which the net charge of the metal ion does not 
change (i.e., contributions to V^oiarize rather than Vhop in 
the language of £111 Ap . one of (A,B) should correspond 
to a Cu state, i.e., metal b\ g state and one to a Pc state, 
e.g., ligand e g state and similarly for (D,E). 

Hence, overall the four indices may involve one entry 
for an e g LUMO state, two entries for the b\ g state, and 
one entry for a doubly-filled ligand state: the highest- 
lying such states are a\ u , a,2 U , or J>2 ligand states for 
single-molecule electronic structure Q. However, be- 
cause a\ u and am are odd under inversion, but b\ g and 
e g are even, the two-electron integrals involving ai u and 
CL2u are zero. Furthermore B^g transforms like xy in D±h 
symmetry, b\ g like x 2 — y 2 , and e gx . y like zx,zy. The 
two-electron integral involving bi g is therefore odd in ei- 
ther y or in x depending which e g state appears. So, in 
fact the only important doubly-occupied states are the 
filled e g states which appear slightly below the a\ u and 

In order to simplify the calculation we assume there 
is only one electron-hole pair produced in the Cu(II)Pc 
dimer (additional electron-hole pairs will cost more en- 
ergy) . As in £1111 A 21 we need consider only the situation 



where M, 



total 



m total 



Ivl total 



in order to extract 



the exchange constant. We find the Hamiltonian v can 
be written as the linear combination of the product of the 



Jl 


Aa 2 t 


(15) 


(U g -U x -E x - 2j eh ) 2 ' 


a 


= 2{X,b\b,G}; 


(16) 




= {G, G\G, G}\ 


(17) 




= {G,X\G,X}; 


(18) 


jeh 


= {G,X\X,G}. 


(19) 



a measures the Cu spin's ability to polarize the ligand, 
U g is the Coulomb interaction between electrons in the 
filled Eg state, U x is the Coulomb interaction between 
electron and hole within one molecule, j h is the electron- 
hole exchange integral, and E x is energy gap between 
LUMO and filled e g state. From equation (fT5|) . we can 
see that the magnitude and sign of J\ depend on the 
inter- molecule transfer integral t. We calculate the ma- 
trix element t for polarization transfer by considering the 
individual hole and electron hoppings among the four 
states below (Figure E]). We find t = j^fc, where t g>e 
are the single-particle transfer integrals between the filled 
states and LUMO of different molecules, respectively. 

If we consider the contributions from both components 
of e g symmetry, the total transfer integral reads 



t x +t v 
u + E x ' 



t x 
t y 
f 

t 



= m 



eg' 

t v t y - 



o A,LUMO 

°g,i 

A, filled 



jt I B,LUMO\. 
n C ore\e. nA ), 



9,i 



o B, filled; 



(20) 

(21) 
(22) 

(23) 

(24) 



where i € {x, y}, H core is the core Hamiltonian for a 
Cu(II)Pc dimer. In the present calculations we evaluate 
H core using the Gaussian 98 code [jjE], using the same 
basis set and exchange-correlation functional described 
above. The symbols A,B refer to these two molecules, 
and ,B ) refer to the e g -symmetry single-molecule lig- 
and states belonging to molecule A or B. 

We use the single-molecule orbital coefficients of the 
isolated molecules and the core Hamiltonian for the 
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FIG. 6: The four states for calculating the transfer integrals 
between polarized triplet states of different molecules as the 
example. 



molecular dimcr to calculate the transfer integrals t e ,t g 
in the molecular configurations with different stacking 
angles (20° — 90°) as shown in Figure [T] The distance 
between these two planes is 3.4 Angstroms [13, EH. In 
FigureEJ we show the variation of inter-molecule hopping 
integrals with stacking angle; t g and t e change both mag- 
nitude and sign with stacking angle; this contributes to 
corresponding changes in the polarization hopping ma- 
trix element t and the exchange constant J-±. 

In Figure [51 we display t x + t y , which contributes the 
dependence on stacking angle to t and hence to J\, as 
a function of stacking angle in the range (20° — 90°). 
When the angle is equal to 45°, we find weak ferromag- 
netic (nearly paramagnetic) coupling. When the angle is 
equal to 65°, the magnetic interaction is relatively strong 
anti-ferromagnetic. This calculation qualitatively agrees 
with the experimental results [l6|, though this calcula- 
tion cannot predict the absolute magnitude of the ex- 
change coupling. 



2. Mn(II)Pc 



3.0 
2.5 
2.0 
1-5 

V 

<D 1.0 

I M 
to 

X 0.0 

>-0.5$- 
■1.0- 
■1.5 
■2.0 



3.5-i 

3.0 

2.5 

2.0 

1.5- 

1.0- 

0.5- 

0.0 
-0.5 
-1.0- 
-1.5- 
-2.0- 



* e 



/V. 



<* 30 40V 50 60 70 80 90 

\ / stacking angle(degrees) 



(a) 



(0 

I 



-* e 



■yN-t ; / , 

^30/ 40 50 60 70/" 80 90 



/ slacking angle(degrees) 



(b) 



FIG. 7: (Color online.) 
defined in equations (f20l 
the filled e g states and (b) t 



The variation in hopping integrals 
(|24|l with stacking angle: (a) t g for 
for the empty e g states. In 



each figure, the solid black curve with square points denotes 
hopping integrals between a;-oriented states of the different 
molecules, and the dashed red curve with triangular points 
shows the integrals between y-oriented states. 



< 




f »T-*~~t- 
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The Mn(II)Pc calculation is more complicated because 
there are three unpaired electrons per molecule which oc- 
cupy a\g, and e g states, so it is necessary to use group 
theory to simplify the calculation of the two-electron inte- 
grals. By a similar procedure to Cu(II)Pc (the details are 
shown in Appendix f5]), we find a weak ferromagnetic in- 



FIG. 8: Variation of the indirect exchange Ji with the stack- 
ing angles shown in Figure[T] arbitrary units (A.U.) are used. 



teraction when the stacking angle is 45° but a relatively 
strong anti- ferromagnetic interaction for 65°. Unfortu- 
nately, even when combined with the super-exchange re- 
sults for Mn(II)Pc obtained in i jlll A 21 (which always pro- 
duce anti- ferromagnetic exchange), this result disagrees 
with the experimental observation of strong ferromag- 
netic coupling near <f> = 45°. 

IV. AB INITIO DFT CALCULATIONS 
A. Cu(II)Pc 

We carry out self-consistent calculations of the elec- 
tronic structure for molecular dimers for the "+" struc- 
tural model at different stacking angles (Figure [p by 
using the Gaussian code with a 6-31G basis set [38| and 
the unrestricted B3LYP (UB3LYP) exchange-correlation 
functional [3^, We perform calculations for differ- 
ent stacking angles ranging from 20° to 90° as shown in 
Figure [Hjl we have tested the convergence of our results 
with respect to basis set by performing a calculation with 
a 6— 31+G* basis set (which includes additional polariza- 
tion functions and diffuse functions) at a single stacking 
angle (45°) and find negligible changes. We compare di- 
rectly the DFT total energies, and hence calculate the 
exchange splitting from the difference of the total ener- 
gies of the broken-symmetry low-spin state and high-spin 
state. For all stacking angles we find it necessary to opti- 
mize carefully the occupancy of the Kohn-Sham orbitals 
in order to ensure that there is no charge disproportiona- 
tion between the molecules; our lowest-energy converged 
states have Mulliken charges of approximately +1.00, and 
nominal spin populations of ±0.68, on each Cu atom. We 
also need to ensure that the numerical convergence error 
in the DFT calculations is much smaller than the order 
of the exchange couplings (IK ~ 10~ 6 Hartree); in our 
calculations, we converge to at least 10~ 9 Hartree. It is 
encouraging that we find negligible spin contamination in 
our final Kohn-Sham wave- functions, i.e. the (S 2 ) com- 
puted for the fictitious non-interacting Kohn-Sham de- 
terminants is close to 2.0 for the triplets ((S 2 ) = 2.0053) 
and to 1.0 for broken-symmetry states ((S 2 ) = 1.0053) — 
note, however^ that this is not the same as the expecta- 
tion value of S 2 in the true many-body wave function. 

In the broken-symmetry state, one b\ g orbital with spin 
up is localized on one molecule; the other with spin down 
on the other molecule as shown in Figure [H Meanwhile, 
in the triplet state, two b\ g orbitals with spin up are 
localized on both molecules. This is consistent with the 
DFT calculation of isolated Cu(II)Pc molecule in which 
localized b\ g state carries the unpaired metal electron. 

The predicted trend of the exchange couplings is con- 
sistent with perturbation theory calculations shown in 
Figure and in particular shows a strong increase in 
the coupling as the molecules approach perfect 7r stack- 
ing (0 = 90°). For the /3 phase, (0 = 45°) we have 
J = E BS - E T = -1.1 x 10" 6 Hartree « -0.3K (See 
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(b) 

FIG. 9: (Color online.) The broken-symmetry orbitals for 
unpaired electrons in each molecule with spin-up (a) and 
spin-down (b) from our DFT calculations for Cu(II)Pc dimer 
(stacking angle 65° ). Notice that 6i 9 states are localized on 
different molecules. 



Figure [T0|) . in agreement with the experimental observa- 
tion of a nearly paramagnetic state at accessible tem- 
peratures, and for a-phase (<f> = 65°) we have J = 
E BS - Br = -5.5 x 10" 6 Hartree w -1.7 K (see Fig- 
ure I10p which gives us a very good agreement with ex- 
perimental observation J ~ — 1.5K. 

The magnitude of the exchange couplings in Cu(II)Pc 
is very small: about 10 -6 Hartree, which is right at the 
edge of the accuracy of the DFT calculation, since there 
will be errors from the imperfect density functionals and 
from the finite basis sets as well as the numerical con- 
vergence errors discussed above. However, we can have 
some confidence in these results for three reasons. First, 
they agree remarkably well with the magnetization mea- 
surements made by the SQUID technique [16J . Second, 
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many of the sources of DFT error could be expected to 
cancel when we compute the energy difference between 
systems that are so similar in every respect except their 
spin orientation. Third, as discussed above, the results 
agree with the trends predicted by perturbation theory. 



■6-1 




1 -I stacking angle(degrees) 



FIG. 10: The energy difference J = Ebs — Et as a func- 
tion of stacking angle from 20° to 90° . Notice the qualitative 
consistency between this figure and Figure [8] 



at 65°, weakly ferromagnetic at 45°), as well as possibly 
other routes. The different mechanisms involve different 
intra-molecular couplings, and so this competition is very 
difficult to quantify on the basis of model calculations. 

Despite the very different methodology, DFT calcula- 
tions on Cu(II)Pc produce results that are remarkably 
consistent with the perturbation theory. When the an- 
gle becomes small, the oscillatory structure of exchange 
interactions calculated by both perturbation theory and 
DFT is a signature of the indirect exchange interaction, 
rather as conventional RKKY oscillations are in a normal 
metal. 
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APPENDIX A: THE INDIRECT EXCHANGE 
FOR MN(II)PC 



V. CONCLUSION AND DISCUSSION 

From perturbative calculations of Cu(II)Pc and 
Mn(II)Pc we find that the exchange interaction between 
two Cu(II)Pc molecules is dominated by indirect ex- 
change. When the stacking angle is 65°, the indirect ex- 
change is predicted to be anti-ferromagnetic, while when 
the stacking angle is 45°, it is very weakly ferromagnetic. 
Both these results agree qualitatively with the experi- 
mental observations (see 

In Mn(II)Pc, by contrast, both super-exchange and in- 
direct exchange contribute. The sign of the indirect ex- 
change interaction in both cases is dependent on the sign 
of inter-molecule electron transfer integrals, and hence 
varies with stacking angle; however, the most important 
terms in the super-exchange are always positive (anti- 
ferromagnetic) . 

The main discrepancy with the experiments is in the 
case of Mn(II)Pc, where our perturbative calculations 
do not give the very strong ferromagnetic interaction 
which was observed. This is probably because the true 
exchange interaction involves the competition between 
super-exchange (always antiferromagnetic) and indirect 
exchange (predicted to be once again anti-ferromagnetic 



First we need to find the symmetry properties of the 
products of pairs of one-electron functions that appear in 
equation ll3l Here we consider the most complicated case, 
the product of two e g states. Eventually, we will consider 
the scattering between filled and empty e g levels in the 
molecule, through interaction with the e g states of the 
Mn ion. To do this, we need the elements of the matrix 
X such that 

X- 1 {e g ®e g ]X = a lg ®b lg ®a 2g ®b 2gi (Al) 

which are the Clebsch-Gordan coefficients for the product 
representation e g ® e g . We can label them as X(a,ij) 
where a refers to one of the irreducible representations 
appearing on the right of equation (| Al|) . and i,j label 
the functions transforming as e g . We find 





Olfl 


a 2g 


b ig 


b 2g 


XX 


(1/V2 





1/V2 





xy 





1/V2 





l/\/2 


yx 





-1/V2 





1/V2 


yy 


UA/2 





-1/V2 






Because [7—771 belongs to the identity representation, we 
can then rewrite v as: 
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-^(M - 1^>) 


(A12) 


= 


- \yy)) 
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-j=(\xy) + \yx)). 


(A14) 



We use X, Y to label the ligand E gXl E gy states, and x, y 
to label Mn e gx , e gy orbitals. Now we introduce operators 
which create electron-hole excitations with different spin 
symmetries on the Pc: 



4 ( °' 0) = at* = -*=(aP+aF); 



t(l,0) _ JT _ 1 



a 



tTJ 

t(i.i) 
t(i.-i) 



= a. 



= — =(a 



tTJ. 



,UT 



^x^sj; 



aJ iT 



); 



,tTT 



a UJ 



c l x T Ci sJ; 
^IxJ^sT' 



(A15) 

(A16) 

(A17) 
(A18) 
(A19) 



where i runs over the two orientations of the ligand e g 
states (i = A, Y) and the subscripts g, x of i label the 
filled ligand e g states and LUMO ligand e g states. The 
following operators characterize the spin degrees of free- 
dom within the subspace where no charge transfer takes 
place: 



S* = 



1, 

°jl c jh 



hi 



(A20) 
(A21) 
(A22) 



where j runs over all the e g states of the Mn ion and the 
ligand: j = X x , X g , Y x , Y g , x, y. Using these operators, 
we can expand v as, 

v = v\ + V2 + i>3 + spin-independent terms, (A23) 
where 

til = (a s x +a^)(2V2(P a + P4-Pi/2-P 2 /2)) 

+2V2a T x {P 1 S z + P 2 S Z ) + 2V2al{P 2 S z x + PiS*) 
+2S|(p ia tU + P 2 < U ) + 2St(p l0 tU + P 2 a^ U ) 

-2S-(Pia^ + P 2 atJ T ) - 2S-(Pi4, TT + p 2 a x V ) (A24) 
V2 = 2{n x] [-P 5 n Xx] - Pencil + n xl [-P5n Xxl - P6«y x |] 
+nj /T [-P 5 ny xT - P 6 n X:cJ ] + n yl [-P 5 n Y:cl -P 6 njf x i] 
+St[-P 5 5 Xx - P 6 Sy J + St[-P 5 Sy x - P 6 S" ] 



+S-[-P5S x -P 6 S. 



-p 5 4 



P 6 5+ ]) (A25) 



D3 



+n yT [-P^ny gT - Penxgtl +%|[-f5 n y g i ~ P 6 n x g |] 



+S|[-P^ X 



ns Xa 



P 6 S Yg 

- PkS v 



- Sy [- P 5 S Y, 



P 6 S Xg] 

- P^SU), (A26) 
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and 



Pi 


= {X,X X \X 9 ,X} 


+ {y,Y x \Y g ,y} 




= {x,Y x \Y g ,x} + 


{y,x x \x g ,y} 


p 3 


= {X X ,x\Xg,x} 


+ {Y x ,y\Y g ,y} 


Pi 


= {X Xl y\X g ,y}- 


f {Y x ,x\Y g ,x} 


p 5 


= {x,X x \X x ,x} + {y,Y x \Y x ,y} 


p 6 


= {x,Y x \Y x ,x}^ 


-{y,X x \X x ,y} 


p' 


= {X,Xg\Xg,X} 


+ {y,Y g \Yg,y} 


P' 


= {x,Yg\Yg, X } + 


{y,X g \X g ,y}. 



To get the leading terms in the effective Hamilto- 
nian [5j we need consider the situation where the to- 
(A27) tal z-direction angular momentum on both molecules is 
M to tai = Mf nA + M s Mns + A/ s Pc = 2. If we restrict 
ourselves to excitations in which there is only only one 
electron-hole pair in total on the two Mn(II)Pc molecules, 
and suppose it resides in the X-symmctry orbitals (the 
Y-symmetry states are completely decoupled in the "+" 
model), we are left with a total of 18 states . These are 
made up as follows: 



Here we can see v\ governs the creation of the electron- 
hole pair, while v 2 and v$ represent the exchange inter- 
actions between spins on the Mn ion and on the ligand. 
Using this form of z), we can build the Hamiltonian matrix 
for two sets of wave functions: those in which the total 
z-component of spin on one molecule (Mn plus ligand) 
is respectively +3/2 and +1/2. We label the individual 
states as |5 , M n (ii): Sp c ), where the first index is the spin 
configuration of Mn ion, and the second is the spin con- 
figuration of the ligand in the X otY spatial component. 

1. The +3/2 states are |(3/2),g), |(3/2),(5 = 0, M = 

0) ), |(3/2), (S = 1,M= 0)), |(l/2), (5 = 1,M = 

1) ), and the corresponding matrix is 








V 



fa 7 






V6S 



7 
y/66 







(A28) 
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where a = 2\/2/3(Fi + P 2 ),(3 
l/2Px-l/2P 2 ),7 = ^/3(P 5 ' + 
l/3(-P^-P^-P 5 -P 6 ). 



2 V / 2(P 3 + Pi 
-P 5 -P 6 ),S 



2. The +1/2 states are |(l/2),g), |(l/2),(5 = 0, M = 

0) ), |(l/2), (S=1,M = 0)>, |(-1/2),(5=1,M = 

1) ), |(3/2), {S = 1,M = -1)), and the matrix of 
Coulomb interaction vi is 




(A29) 





= 3/2,M s Mns = 


1/2, Mf c = 


0)- 


5 states ( 


(M™ nA 


= l/2,M a MllB = 


3/2, M s Pc = 


0)-» 


5 states 


(M™ nA 


= l/2,M s Mns = 


1/2, A/ a Pc = 


1)- 


2 states 


(M™ nA 


= 3/2,A/ a Mns = 


-1/2, Mj c 


= 0) 


— » 2 states 


(M™ nA 


= -l/2,M s MnB 


= 3/2, Mf c 


= 0) 


— » 2 states 


(M™ nA 


= 3/2,M a Mns = 


3/2, Aff c = 


-1) 


— » 2 states. 



The perturbation includes both the Coulomb interaction 
v discussed above, and the hopping t which transfers an 
electron-hole pair from one Mn(II)Pc to another. We 
found the leading term of spin-| couplings J\ to be: 



Ji 



3a 2 t 



{U g -U X ~E X - 2 Joh ) 2 



(A31) 



where the definitions of U g , U. 



t P and t g are the 



Px i jch , 

same as those in the Cu(II)Pc calculation 

Considering other situations such as the product of 
aig and e g states gives qualitatively similar results that 
depend in the same way on the transfer integrals between 
molecules. Finally, we include excitations through both 
components (X and Y) of the ligand e g states, so as in 
the Cu(II)Pc calculation we should combine the electron- 
hole pair transfer integrals to form t = t x + t y . 
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